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Abstract 

The  virtual  test  bed  (YTB)  computation  environment  provides  an  easy-to-use  way  for  electric  system  modeling  and  simulation.  Resistive 
companion  (RC)  modeling  of  batteries  in  YTB  is  presented  in  the  paper.  Native  RC  battery  modeling  approach  and  the  one  using  current 
controlled  voltage  source  (CCVS)  interface  are  presented  and  compared  through  two  battery  models  with  different  degrees  of  complexity. 
The  RC  battery  models  are  validated  by  comparing  VTB  simulation  data  to  those  calculated  directly  through  stand-alone  Fortran  codes.  It  is 
concluded  that  using  CCVS  interface  is  currently  the  first  choice  in  modeling  batteries  in  VTB.  Simulations  using  the  RC  battery  models  in 
VTB  are  also  presented  and  analyzed.  It  is  shown  that  RC  modeling  provides  a  powerful  way  for  the  simulation  of  battery  systems  in  VTB. 
©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

In  the  last  decade,  the  development  of  electric  and  hybrid 
vehicles  and  the  rapid  expansion  of  portable  computers,  con¬ 
sumer  electronics  and  telecommunication  tools  in  the  world 
have  spurred  research  into  advanced  battery  systems  with  em¬ 
phasis  on  operation  under  various  conditions  and  on  materials 
and  cell  construction  to  enhance  cycle-life  and  performance. 
The  research  effort  has  been  undertaken  by  both  extensive 
experimental  tests  and  computer  simulations  based  on  mod¬ 
els  evolved  from  physical  and  chemical  laws  of  the  processes 
occurring  in  the  cell,  i.e.,  first  principles  [1]  or  empirical  mod¬ 
els.  These  models  play  an  important  role  in  cell  design  in  pre¬ 
dicting  cell  performance  under  various  conditions,  and  more 
recently,  in  the  integration  into  electric  system  models  [2,3] 
to  highlight  the  real-time  behaviors  of  the  whole  system.  The 
latter  feature  is  usually  accomplished  with  circuit  simulators, 
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such  as  SPICE  and  virtual  test  bed  (VTB)  [2,4,5]  developed 
at  University  of  South  Carolina. 

The  VTB  software  is  dedicated  to  provide  a  unique  capa¬ 
bility  for  simulation  and  virtual  prototyping  of  power  elec¬ 
tronic  systems.  The  major  modules  of  VTB  software  are 
schematic  editor,  visualization  engine,  system  solvers  and 
entities.  VTB  entities  are  the  models  of  particular  devices, 
i.e.,  resistor,  capacitor  and  batteries,  etc.  They  are  coded  in 
resistive  companion  (RC)  [4-7]  format  which,  by  discretiz¬ 
ing  and  linearizing  time  differential  equations  into  algebraic 
equations,  represents  the  physical  circuit  through  a  set  of  si¬ 
multaneous  algebraic  equations  analogous  to  those  describe 
a  direct  current  (dc)  circuit.  In  the  RC  method,  entities  are 
connected  to  each  other  through  ports  or  terminals.  Each 
RC  entity  only  handles  its  own  governing  equations.  VTB 
system  solvers  collect  the  feedback  from  entities,  that  is, 
conductance  matrices  and  current  source  vectors  to  generate 
system  matrices  and  solve  all  interconnected  devices  in  the 
system. 

Batteries  and  fuel  cells  are  important  power  sources  in 
the  electric  system  simulation.  In  the  following,  the  paper 
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Nomenclature 

a  specific  interfacial  area  in  the  cathode  of  the 

Li-ion  battery  (m2  m-3) 

ane g  specific  surface  area  of  negative  electrode  of 

the  NiMH  battery  (cm2  cm-3) 
apos  specific  surface  area  of  positive  electrode  of  the 
NiMH  battery  (cm2  cm-3) 

Aneg  geometry  area  of  negative  electrode  of  the 
NiMH  battery  (cm2  cm-3) 

Apos  geometry  area  of  positive  electrode  of  the 
NiMH  battery  (cm2  cm-3) 
b  RC  current  source  vector  element 

B  RC  current  source  vector 

c  concentration  of  electrolyte  in  the  Li-ion  bat¬ 

tery  (mol  m  3) 

ce  concentration  of  KOH  electrolyte  in  the  NiMH 
battery  (mol  cm-3) 

cmax  maximum  Li  salt  concentration  in  polymer  in 
the  Li-ion  battery  (molm-3) 
cmh  concentration  of  hydrogen  in  metal  hydride  in 
the  NiMH  battery  (mol  cm-3) 

CNi(OH)o  concentration  of  Ni(OH)2  in  the  NiMH  battery 
(mol  cm-3) 

cs  concentration  of  lithium  ions  in  the  solid  phase 

of  the  Li-ion  battery  (molm-3) 
cs, average  average  lithium  ion  concentration  in  an 
intercalation  particle  in  the  Li-ion  battery 
(molm-3) 

cj  maximum  concentration  of  lithium  ions  in 
solid  in  the  Li-ion  battery  (molm-3) 

D  diffusion  coefficient  of  electrolyte  in  the  poly¬ 
mer  in  the  Li-ion  battery  (m2  s-1) 

Ds  diffusion  coefficient  of  lithium  ion  in  the  solid 
matrix  in  the  Li-ion  battery  (m2  s-1) 

Eq  open-circuit  potential  of  the  battery  (V) 
f±  mean  activity  coefficient  of  the  electrolyte  in 
the  Li-ion  battery 

F  faradic  constant  (98647  C  equiv. - 1 ) 

g  RC  conductance  matrix  element 

G  RC  conductance  matrix 

h  time  step  (s) 

i'oi  exchange  current  density  of  anode  reaction  in 
the  Li-ion  battery  (Am-2) 

iok  exchange  current  density  of  reaction  k  in  the 
NiMH  battery  (A  cm-2) 

/cell  applied  current  (A) 

4  current  at  node  k  in  the  RC  model  (A) 

I  superficial  current  density  in  the  Li-ion  battery 

(A  m-2) 

jn  pore  wall  flux  of  lithium  ions  (mol  m-2  s) 

&2  reaction  rate  constant  at  cathode/polymer  in¬ 
terface  in  the  Li-ion  battery  (m4  mol-1  s-1) 


/c  thickness  of  the  cathode  in  the  Li-ion  battery 

(cm) 

/neg  thickness  of  negative  electrode  in  the  NiMH 
battery  (cm) 

/pOS  thickness  of  positive  electrode  in  the  NiMH 
battery  (cm) 

/s  thickness  of  the  separator  in  the  Li-ion  battery 

(cm) 

A;, neg  equivalent  thickness  of  metal  hydride  material 

in  the  NiMH  battery  (cm) 

/^pos  equivalent  thickness  of  nickel  active  material 
in  the  NiMH  battery  (cm) 

Lmh  loading  of  metal  hydride  material  in  the  NiMH 

battery  (gem-2) 

Lni(oh)2  loading  of  nickel  active  material  in  the  NiMH 
battery  (gem-2) 

Mmh  molecular  weight  of  metal  hydride  mate¬ 
rial^  mol-1) 

MNi(OH)o  molecular  weight  of  nickel  active  material 
(gmol-1) 

n  number  of  electrons  transferred  in  reactions  in 

the  Li-ion  battery 

Ni  molar  flux  in  v-direction  of  species  i  in  the  Li- 
ion  battery  (mol  m-2  s) 

p  gas  pressure  in  the  NiMH  battery  (atm) 

r  distance  normal  to  surface  of  cathode  material 

in  the  Li-ion  battery  (m) 

R  ideal  gas  constant  (8.3145  J  mol-1  K-1) 

RCell  battery  resistance  (Q) 

Rs  radius  of  cathode  material  in  the  Li-ion  battery 

(m) 

s+  stoichiometric  coefficient  of  Li  ions  in  the  Li- 
ion  battery 

SOC  state  of  charge  in  the  NiMH  battery 
t  independent  time  variable  (s) 

transference  number  of  lithium  ions  in  the  Li- 
ion  battery 

T  temperature  (K) 

f/eq,/,ref  equilibrium  potential  of  reactions  /  at  reference 
reactant  concentration  in  the  NiMH  battery  (V) 
vceii  close-circuit  potential  of  the  battery  (V) 

Vk  node  potential  at  node  k  in  the  RC  model  (V) 

v+  number  of  anions  into  which  a  mole  of  elec¬ 

trolyte  dissociates  in  the  Li-ion  battery 
Vgas  gas  volume  in  the  NiMH  battery  (cm3) 

Lout  output  potential  obtained  in  battery  SFM  model 

(V) 

v  spatial  coordinate  in  the  Li-ion  battery  (m) 

z+  charge  number  of  lithium  ions  in  the  Li-ion 
battery 

Greek  letters 

A0neg  potential  difference  at  the  solid  and  liquid  in¬ 
terface  on  the  negative  electrode  in  the  NiMH 
battery  (V) 
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A0pos 

potential  difference  at  the  solid  and  liquid  in¬ 
terface  on  the  positive  electrode  in  the  NiMH 
battery  (V) 

£ 

porosity  of  electrode  in  the  Li-ion  battery 

0 

cathode  overpotential  in  the  Li-ion  battery  (V) 

K 

conductivity  of  electrolyte  in  the  Li-ion  battery 
(Sm-1) 

PMH 

density  of  metal  hydride  active  material  in  the 
NiMH  battery  (g  cm-3) 

PNi(OH)2 

density  of  nickel  active  material  in  the  NiMH 
battery  (gem-3) 

O 

conductivity  of  sold  matrix  in  the  Li-ion  battery 
(S  m-1) 

<t>\ 

solid  phase  potential  in  the  Li-ion  battery  (V) 

02 

liquid  phase  potential  in  the  Li-ion  battery  (V) 

Superscript 

0 

initial  condition 

presents  and  compares  RC  battery  modeling  methods  which 
evolve  with  the  complexity  of  battery  models. 


2.  Resistive  companion  modeling  of  batteries 


Fig.  1.  Battery  RC  model  with  Thevenin-Norton  transformation. 


where  the  conductance  matrix  is  of  the  form  G(t  —  h)  = 

goo  goi 


g  10  gll 

bo 


,  the  current  source  vector  of  the  form  B(t  — 


h)  = 


b\ 


,  t  the  time  and  h  the  time  step.  The  conductance 

matrix  elements  are  given  by 

di 


gjk(t-h)  = 


j 


where  j,  k  =  0,  1 


dvk)t-h 

and  current  source  vector  elements  are  given  by 
bj(t  -  h)  =  —ij(t  —  h)  +  Y  gjkVkif  -  h),  where 


(5) 


j,  k  =  0,  1 


(6) 


2.7.  Overview 

A  simple  battery  model  as  described  by  Eq.  (1)  is  referred 
to  explain  RC  battery  modeling: 

Fcell  -  f  0  f  Cell^cell  (1) 

where  vceii  is  the  close-circuit  potential,  Eo  the  open-circuit 
potential,  /ceii  the  applied  current,  and  Rce n  the  internal  resis¬ 
tance  of  the  battery.  In  Eq.  (1),  it  is  assumed  that  the  battery 
has  a  constant  open-circuit  potential  and  close-circuit  poten¬ 
tial  varies  linearly  with  applied  current.  Such  a  battery  is  often 
viewed  as  a  voltage  source  with  internal  resistance.  However, 
resistive  companion  method  requires  all  entities  to  be  repre¬ 
sented  as  current  source  and  resistance  parallel  combinations 

[5] .  The  battery  needs  to  be  replaced  with  an  equivalent  cir¬ 
cuit  consisting  of  the  parallel  combination  of  a  resistor  and 
current  source  by  applying  Thevenin-Norton  transformation 

[6]  (see  Fig.  1).  The  close  circuit  voltage  vceii  and  applied 
current  /Cell  can  be  related  to  terminal  potential  and  current 
variables  by 

Fcell  =  Vo  -V\  (2) 

Ae  ll  =  io  =  ~h  (3) 

where  vo,  v\  are  node  potentials  and  io,  i\  the  node  currents. 

For  an  electric  device  with  two  terminals,  the  resistive 
companion  method  expects  the  following  relation: 

(4) 


From  Eqs.  (l)-(6),  we  can  obtain  the  RC  representation 
for  the  battery: 


*oO) 

r  l 

^cell 

l  i 

^cell 

vo  (t) 

r  Eo  i 

PceW 

i  i  (t) 

1 

1 

Vi  (t) 

Eo 

-  ^cell 

^cell  - 

-  Pcell  - 

The  conductance  matrix  and  current  source  vector  are  used 
by  VTB  solvers  in  the  numerical  calculations.  If  the  battery 
has  fixed  open-circuit  potential  Eo  and  internal  resistance 
Rceii  which  is  the  case  in  Eq.  (1),  the  elements  of  its  con¬ 
ductance  matrix  and  current  source  vector  are  all  constants. 
Devices  of  such  kind  are  called  linear  entities  [5]  in  VTB. 
However,  the  battery  model  (Eq.  (1))  is  oversimplified,  i.e.  it 
does  not  reflect  the  change  of  the  open-circuit  potential  with 
the  state  of  charge  (SOC)  and  the  battery  internal  resistance 
may  also  vary  in  the  process  of  discharge.  Many  battery  mod¬ 
els  available  in  the  literature  are  more  sophisticated  than  Eq. 
(1)  and  are  usually  represented  by  a  set  of  nonlinear  algebraic- 
differential  equations  yielding  nonlinear  entities  [5]  in  VTB 
environment.  In  the  following  sections,  two  battery  models 
with  different  complexity  chosen  from  the  literature  are  dis¬ 
cussed  in  details  as  extensions  of  the  RC  battery  modeling  in 
VTB  environment. 

2.2.  Native  RC  modeling  of  a  nickel  metal  hydride 
(NiMH)  cell 

A  nickel  metal  hydride  cell  model  described  by  Wu  et  al. 

[7]  is  used  here  to  introduce  native  RC  modeling  of  batteries 


7(0  =  G(t  -  h)V(t)  -  B(t  -  h ) 
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of  which  the  open-circuit  potential  varies  with  the  SOC.  It 
is  assumed  in  the  model  that  both  the  positive  and  negative 
electrodes  are  planar  electrodes  and  diffusion  process  inside 
the  solid  active  material  is  not  important.  Four  reactions  in 
the  nickel  metal  hydride  cell  have  been  considered.  They  are: 

•  Nickel  electrode 


discharge 

NiOOH  +  H20  +  e"  <=*  Ni(OH)2  +  OH 

charge 


1 


-02  +  H20  +  2e 


20H“  - 
Metal  hydride  electrode 


discharge 

MH  +  OH-  <=►  H20  +  M  +  e 

charge 


-02  +  H20  +  2e 


20H 


/ 


The  governing  equations  for  the  cell  are: 
dCNi(OH)2  j  1 


/ 


y,pos  d, 

dCMH 


F 


73 

F 


y’neg  dt  - 
f^gas  d/?o2  ^pos^pos^pos72  T  r*neg4iegTlneg74 

hr  d  t  f 

*cell  —  ^pos^pos^pos(jl  T  72) 

—  *cell  =  ^neg^neg^neg(73  T  7*4) 
t’cell  =  (^0pos  —  A0neg)  T  *cell^cell 


where 


7i  =  «oi 


(  cNi(OH)2 
\CNi(OH)2,ref. 


Cf 


ce,ref 


x  exp 


0.5  F 
~RT 


(A0pos  F^eq,l,ref) 


-  1- 


CNi(OH)2 
CNi(OH)2,ref 

0.5  F 


x  exp 


RT 


(A0pos  F^eq,  l,ref) 


72  =  *02 


(  Af )  exP  (  RT(Mpos  £4q, 2, ref) 


PO2 
/-’(Kiel 


0.5 


F 


X  exp  |  (A(ppOS  f/eq,2,ref) 


(8) 


(9) 


(10) 


(11) 


(12) 

(13) 

(14) 

(15) 

(16) 
(17) 


(18) 


(19) 


P  o2  . 

J4  —  *4, ref 

PO2 ,  ref 

,  _  LNi(OH)2 

ly, pos  —  j 

PNi(OH)2  *pos*2pos 

j  _  Lmh 

hi  neg  — 

PMH*pos*2pos 

The  initial  conditions  for  Eqs.  (12)— (14)  are: 


(21) 

(22) 

(23) 


Ni(OH)2  =1_SQC0 
CNi(OH)2,max 

CMH  =  SQC0  ^Ni(QH)2ApOSMMH 
CMH.max  ^MH  AnegAfNi(OH)2 

Pq  =0.1  atm 


(24) 

(25) 

(26) 


Eqs.  (12)— (17)  need  to  be  converted  to  RC  format  by 
discretizing  equations  within  one  time  step  and  calculating 
the  conductance  matrix  and  current  source  vector  using  Eqs. 
(5)  and  (6).  As  seen  in  Eqs.  (12)— (17),  the  model  equations 
are  nonlinear  with  respect  to  the  variables  such  as  A0pos, 
Qsfi(OH)2  pmh?  and  A0neg.  Analytical  solutions  for  the  RC 
matrices  are  not  available.  One  common  practice  is  to  use 
finite  difference  approximation  to  calculate  the  RC  matrices 
numerically: 


Sjk  ~ 


=  (_  |  yi+k  d<ce" 

t-h  dnCell 


,  .  y+k  *cell  T  A i  f cell 


where  j,  k  =  0,  1 


(27) 


,  (  dz 

bj  =  -h  +  zlsjkVk  =  (-1  y  t- 


cell 


\  dnceii 


^cell  *cell 


.  \j  (  *cell  T  A  i  *Cell 


=  (-iy 


V  ^cell  T  An  Uceii 
where  j,  k  =  0,  1 


^cell  *cell 


(28) 


The  perturbation  of  the  cell  voltage  (Av)  or  cell  current 
(A/)  must  be  small  enough  to  provide  good  numerical  accu¬ 
racy  for  the  RC  matrices.  A  Fortran  solver  call  GNES  is  used 
to  provide  model  solutions  numerically. 

2.3.  RC  modeling  of  a  Li-ion  cell 

First  principle  battery  models  are  built  on  our  under¬ 
standing  of  the  battery  physics,  including  the  microscopic 


73  =  *03 


CMH 


C( 


CVIH,ref  /  \^e,ref 

0.5E 

exp 


exp 


0.5E 

-RT 


(A0neg  f^eq,3,ref) 


RT 


(A0n  £7eq?3?ref) 


(20) 
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current  composite  polymer  lithium  foil 

collector  cathode  electrolyte  anode 


90?-<2>i)_/  a  80 1  RT  (  s+  t\  \  dine 

dx  k  k  dx  F  \nv+  Z+v+  I  dx 

(35) 


s+  d2@\ 
nF  dx 2 


(36) 


jn  =  k2(cmax  -  C)°-5C°'5 

Cs\r=Rs  exp  U')\ 

X  ^  ' 

_-(cT  -  Cs|r=^s)exp  (-^(>7  -  U')) 


(37) 


Fig.  2.  Lithium/polymer  cell  sandwich. 

electrode  reactions  and  transport  phenomena  inside  the 
porous  electrodes.  Some  first  principle  models  further  in¬ 
corporate  coupled  electrochemical  and  thermal  behaviors  to 
reflect  the  thermal  influence  on  the  performance  of  the  bat¬ 
tery.  They  are  more  rigorous  in  describing  battery  behaviors, 
but  are  more  complicated.  A  first  principle  battery  model 
described  by  Doyle  et  al.  [8]  is  presented  to  illustrate  how 
to  model  such  complicated  batteries  through  RC  method  in 
VTB  environment. 

Doyle  et  al.  [8]  used  the  model  to  simulate  the  galvano- 
static  charge/discharge  behavior  of  the  cell  sandwich  shown 
in  Fig.  2.  The  model  considers  one-dimensional  transport  of 
lithium  ion  from  the  lithium  anode  through  the  polymer  sep¬ 
arator  into  the  composite  cathode.  The  reaction  at  the  anode 
is  assumed  to  take  the  form: 

Li  +  (9p  ^  Li+  —  (9p  +  e-  (29) 

And  the  reaction  at  the  cathode  is  assumed  to  take  the 
form: 


dc  s 
~dt 


=  A 


r  d2cs  2  3 cs 

dr 2  r  dr 


with  the  boundary  conditions: 


dc 


,0 


~Ddl+It+jlx=k+h-  F 


dc 

dx 

dCs 

dr 


x 


=o  =  0 


lr=0  =  0 


9cs 

jn  —  —F)s——\r=R 
dr 

30i  I 

—  U=o  = 

dx  a 


30i 

dx 


I  x=lc—  —  6 


and  the  battery  close  circuit  voltage  is  defined  as 


(38) 

(39) 

(40) 

(41) 

(42) 

(43) 

(44) 


Li+  -  <9p  +  (9S  +  e  <=*  Li+  -  <9S  +  <9p 


(30)  nceii  —  0i|x=o  0iU=/c+zs 


(45) 


The  governing  equations  for  the  Li-ion  battery  are: 
at  the  anode: 

0.50(01 1*=/  +/  -  <p2\x=lc+h) 


I  =  i  01 


exp 


exp 


RT 

0.50(01  L=/c+/s  -  ^2|jc=Zc+Zs) 


RT 


in  the  separator: 

dc  3  2c 
—  =  D 
3 1 


3 1% 


dx2  z+v+F  dx 


302  kRT  I 

I  =  — K  — —  —  -  ( - 1- 


r° 


31nc 


dx 


0  \  nv+  Z+v+  /  dx 


in  the  composite  cathode: 


dc  ic  dzc  I  +  <j(30i/3v)  3z5.  ajn(l  —  tF) 
s—  ~  s  D — - - —  H - — 


(31) 


(32) 


(33) 


3 1 


dx2 


z+v+F 


dx 


(34) 


The  equations  mentioned  above  to  model  the  Li-ion  cell 
is  computation  intensive  because  the  solutions  for  both  x- 
direction  and  r-direction  need  to  be  obtained.  Venkat’s  ap¬ 
proximation  [9]  is  applied  to  eliminate  the  need  to  solve  Eq. 
(38)  in  r-direction  by  converting  the  equation  to  two  differ¬ 
ential  and  algebraic  equations: 


average  ,  ,  \  r\ 

H  (CSj  average  C  s\r=Rs)  —  tl 


3 1 


R2 

s 


5DS 

Rs 


(c 


s,  average 


Cs 


(46) 

(47) 


To  solve  the  model  equations,  finite  difference  approxima¬ 
tion  is  used  to  discretize  the  special  coordinate  x.  Hundred 
internal  node  points  are  used  to  achieve  good  convergence 
and  desired  accuracy.  After  discretization,  the  total  number 
of  equations  is  close  to  103.  A  Standard  Fortran  solver  called 
DASSL  is  used  to  provide  numerical  solution.  If  the  Li-ion 
battery  is  modeled  in  the  native  RC  method,  the  RC  model  is 
a  nonlinear  entity  in  VTB  and  must  be  called  several  times 
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in  each  time  step  in  order  for  system  to  converge  upon  a  so¬ 
lution  at  time  t  [5].  And  in  each  call,  model  equations  have 
to  be  solved  twice  to  get  battery  resistance.  Therefore,  the 
native  RC  modeling  method  is  not  expected  to  be  efficient 
in  dealing  with  such  batteries  where  the  number  of  battery 
model  equations  is  large,  i.e.,  103.  It  is  desired  that  a  new  RC 
modeling  method  could  be  used  to  achieve  higher  numerical 
efficiency  than  the  native  one. 

It  should  be  noticed  that  the  close  circuit  voltage  ob¬ 
tained  through  Eq.  (45)  in  the  Li-ion  battery  model  is  ac¬ 
tually  decoupled  from  the  rest  of  governing  equations  if  the 
charge/discharge  current  is  selected  as  the  control  variable  in 
solving  the  discretized  model  equations.  In  such  case,  there 
exists  a  sparse  Jacobian  matrix  and  numerical  efficiency  can 
be  greatly  improved  by  taking  advantage  of  the  banded  struc¬ 
ture  of  the  Jacobian  matrix  [10].  However,  if  the  voltage  is 
selected  as  the  control  variable,  Eq.  (45)  must  be  coupled  to 
the  rest  of  model  equation  to  solve  the  model  and  the  com¬ 
putation  time  will  be  great  due  to  the  un-banded  structure  of 
the  Jacobian  matrix.  Hence  it  would  be  much  better  that  the 
new  modeling  method  could  exploit  the  decoupled  feature  of 
the  cell  voltage  equation  in  solving  the  battery  model. 

Such  a  RC  modeling  method  using  current  controlled  volt¬ 
age  source  (CCVS)  interface  is  brought  forward  by  the  De¬ 
partment  of  Electrical  Engineering  at  the  University  of  South 
Carolina.  The  method  was  originally  used  to  interface  VTB 
simulator  with  the  ACSL  (advanced  circuit  simulation  lan¬ 
guage)  based  models  [11].  Before  it  is  introduced,  some  re¬ 
lated  background  information  is  provided. 

VTB  simulator  uses  resistive  companion  method  while 
other  simulators  such  as  ACSL  use  signal  flow  method 
(SFM).  When  the  battery  model  is  solved  with  a  Fortran 
solver  like  DASSL,  it  acts  like  a  SFM  model.  For  example, 
with  some  input  of  applied  current  (/in),  the  battery  model  can 
calculate  close  circuit  potential  (Vout)  and  other  state  vari¬ 
ables  reflecting  battery  conditions  as  output  (see  Fig.  3).  An 
interface  or  wrapper  has  to  be  used  if  one  wants  to  connect 
VTB  RC  model  and  SFM  models  together. 

The  CCVS  interface  [11]  used  to  connect  ACSL  and  VTB 
is  adapted  in  the  new  RC  battery  modeling  method.  The  in¬ 
terface  implements  natural  coupling  which  requires  enforce¬ 
ment  of  physical  conservation  principles  to  connect  the  SFM 
battery  model  with  the  RC  battery  model  (VTB  entity)  (see 
Fig.  4). 

The  RC  format  for  the  CCVS  interface  is  formulated  using 
modified  node  analysis  [6] : 
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Fig.  3.  Schematic  of  SFM  type  battery  models. 


Battery  RC  model 


Fig.  4.  Battery  RC  model  with  CCVS  interface,  /is  branch  current  and  I-m  =  /. 


Comparing  to  the  native  modeling  method,  a  branch  cur¬ 
rent  I  is  introduced  as  an  additional  variable  and  the  conduc- 

0  0  1 


tance  matrix  G  = 


0  0  -1 

1  -1  0 


and  current  source  vector 


B  = 


0 

0 

Fout 


are  augmented  to  reflect  the  change.  The  algo¬ 


rithm  works  in  the  following  way.  The  RC  battery  model 
obtains  the  known  values  of  node  voltage  vo,  v\  and  branch 
current  I  at  time  step  t'  -t  from  the  system  solver  and  passes 
the  branch  current  to  the  SFM  battery  model  which  returns 
the  close  circuit  voltage  Vout  as  output.  Then  the  RC  battery 
model  assembles  the  RC  matrices  and  returns  them  to  the  RC 
system  solver  to  determine  new  values  of  node  voltage  vo,  v\ 
and  branch  current  I  at  time  step  f  =  t  +  h. 

Connected  through  the  interface,  the  battery  SFM  model 
appears  as  a  current  controlled  voltage  source  (CCVS)  which 
is  a  linear  entity  in  the  VTB  system.  Consequently,  the  sys¬ 
tem  matrices  are  not  susceptible  to  any  nonlinearity  wrapped 
inside  the  SFM  models  and  modelers  can  select  appropriate 
solvers  for  the  SFM  models  without  affecting  VTB  system 
solvers.  And  most  importantly,  the  battery  SFM  model  needs 
to  be  solved  only  once  in  each  time  step  due  to  the  feature  of 
linear  VTB  entity,  which  is  critical  in  speeding  up  the  simu¬ 
lation. 


3.  Model  validation  and  simulations 

The  VTB  RC  models  are  required  to  be  coded  in  C  lan¬ 
guage  and  implemented  as  dynamic  link  libraries  (DLLs).  A 
DLL  is  an  executable  file  that  acts  as  a  shared  library  of  func¬ 
tions  with  imports  and  exports.  But  it  is  more  efficient  to  code 
the  battery  SFM  models  in  Fortran  than  in  C  because  of  the 
huge  body  of  robust  Fortran  numerical  algorithms  available 
in  public  domain,  such  as  the  standard  DASSL  subroutine 
for  differential-algebraic  equations.  The  battery  SFM  models 
used  in  this  paper  are  also  implemented  as  DLLs.  Such  mod¬ 
ular  structure  can  be  easily  debugged  separately  and  adapted 
to  other  battery  models.  In  order  for  battery  SFM  model  (For¬ 
tran  DLL)  to  correctly  communicate  with  the  CCVS  interface 
and  RC  model  (C  DLL),  mixed  language  programming  [12] 
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Fig.  5.  Simulation  schematic  used  in  model  validation. 

is  required.  The  calling  convention  of  SFM  model  (Fortran 
DLL)  has  to  match  that  of  RC  model  (C  DLL)  and  the  argu¬ 
ments  type  of  Fortran  subroutine  must  be  correctly  declared 
in  the  caller  RC  model.  Both  DLLs  must  be  provided  to  VTB 
to  make  battery  models  work  in  VTB  environment. 

The  validity  of  the  RC  battery  models  is  evaluated  using 
following  simulation  schematic  (see  Fig.  5).  The  PL  in  the 
figure  is  a  programmable  load.  Two  options  are  used  in  the 
validation:  battery  discharging  with  constant  current  load  and 
resistance  load.  RC  NiMH  battery  models  and  RC  Li-ion 
battery  models  using  both  native  modeling  method  and  the 
one  with  CCVS  interface  are  validated  in  the  simulation.  The 
parameters  for  the  battery  models  are  listed  in  Tables  1  and  2, 

Table  1 


Parameters  for  the  NiMH  battery  model 


Parameter 

Value 

flpos  (cm2  cm-3) 

4000.0 

flneg  (cm2  cm"3) 

1000.0 

Apos  (cm2) 

325.0 

Aneg  (cm2) 

360.0 

ce  (mol  cm-3) 

7.0  x  10"3 

Ce,ref  (mol  cm-3) 

1.0  x  10“3 

C'MH.max  (molcm-3) 

1.0  x  10“' 

CMH,ref  (mol  cm"3) 

0.5 

CNi(OH)2,max  (mol  Cm  3) 

3.7  x  10-2 

CNi(OH)2,ref  (molcm-3) 

0.5 

*01  (A  cm-2) 

1.0  X  10-4 

*02  (A  cm-2) 

1.0  x  10"10 

*03  (A  cm-2) 

3.8  x  10“4 

io4  (A  cm-2) 

1.0  x  10"4 

4>os  (cm) 

3.3  x  10"2 

4ieg  (cm) 

2.8  x  10"2 

CNi(OH)2  (g  cm-2) 

6.8  x  10“2 

LMh  (gem-2) 

1.13  x  10“‘ 

AfNi(OH)2  (gmol-1) 

92.71 

Mmh  (g  mol-1) 

70.58 

P02,ref  (<ltm) 

1.0 

Rcell  (^2) 

1.0  X  10“3 

soc° 

0.95 

T  (K) 

298.15 

Ceq,l,ref  00 

0.527 

Ceq,2,ref  00 

0.4011 

Ceq,3,ref  (V) 

-0.8279 

Lgas  (cm3) 

1.01 

PNi(OH)2  (gem-3) 

3.4 

Pmh  (gern-3) 

7.49 

Table  2 


Parameters  for  the  Li-ion  battery  model 


Parameter 

Value 

a  (m2) 

1.0 

Cmax  (molm-3) 

3920.0 

cj  (molm-3) 

29000.0 

D  (m2  s-1) 

7.5  x  10-12 

Ds  (m2s-‘) 

5.0  x  10“13 

*01  (Am-2) 

12.6 

k2  (m4  s- 1  mol) 

1.0  X  10“10 

lc  (m) 

100.0  x  10“6 

4  (m) 

50.0  x  10“6 

n 

1.0 

Rs  (m) 

1.0  X  10“6 

T  (K) 

373.2 

v+ 

1.0 

£ 

0.3 

cr  (S  m-1) 

1.0  X  104 

k  (S  m-1) 

2.176 

respectively.  The  VTB  simulation  results  are  compared  with 
those  obtained  by  solving  system  (Fig.  5)  integratedly  with 
stand-alone  Fortran  codes. 

Simulation  results  are  presented  in  Figs.  6-9.  The  differ¬ 
ences  among  the  results  obtained  through  three  methods  are 
nearly  indistinguishable  as  the  plots  show,  where  the  lines 
are  overlapped.  However,  there  is  a  great  discrepancy  in  the 
computation  time  used  in  the  simulation  as  shown  in  Table  3. 

It  can  be  seen  from  Table  3  that  solving  the  system  with 
stand-alone  Fortran  codes  is  faster  than  using  either  RC  mod¬ 
els.  But  such  an  advantage  is  outweighed  by  the  disadvan¬ 
tage  of  the  need  to  re-program  the  Fortran  codes  when  the 
simulated  system  changes.  It  is  impossible  for  a  stand-alone 
Fortran  code  to  consider  all  possible  combinations  of  large 
variant  models  which  may  be  used  in  the  simulation.  Model¬ 
ing  and  simulation  in  VTB  allows  for  reuse  of  existing  codes 
(entities)  and  complex  system  simulation  can  be  more  easily 
conducted.  It  is  apparent  that  battery  RC  models  in  CCVS 
interface  can  run  more  efficiently  and  obtain  fairly  accurate 


Fig.  6.  Simulation  run  1:  NiMH  battery  discharge  with  constant  resistance 
load  ofO.lf2.Slis  obtained  using  native  RC  battery  model,  S2,  using  CCVS 
interface  model  and  S3  using  stand-alone  Fortran  code. 
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Table  3 


Highlights  of  simulation  results 


NiMH  battery  model 

Li-ion  battery  model 

SI 

S2a 

S3 

SI 

S2a 

S3 

Percentage  errorb  (Vceii) 

Run  1 

<0.005 

<0.01 

N/A 

<0.1 

<0.1 

N/A 

Run  2 

<0.005 

<0.01 

N/A 

<0.01 

<0.05 

N/A 

Simulation  time0  (s) 

Run  1 

~1.2 

—0.4 

-0.2 

>100 

-1.2 

Run  2 

-1.2 

—0.4 

-0.2 

>100 

—0.4 

-0.2 

a  Data  which  have  not  converged  are  excluded  from  calculating  the  percentage  error. 

b  The  percentage  errors  are  calculated  by  comparing  the  cell  voltage  data  simulated  by  the  RC  models  to  the  data  simulated  by  the  stand-alone  Fortran  codes. 
c  Simulation  time  is  based  on  1000  step. 


Fig.  7.  Simulation  run  2:  NiMH  battery  discharge  with  constant  current  load 
of  1  A.  See  Fig.  6  for  the  meaning  of  the  legend. 

results  as  compared  with  those  modeled  in  native  RC  method. 
The  computation  time  saved  from  the  linear  feature  of  RC 
model  with  CCVS  interface  is  significant  when  the  number 
of  model  equations  is  very  large.  Therefore,  it  is  concluded 
from  our  work  that  using  the  CC V S  interface  currently  should 
be  the  first  choice  in  RC  battery  modeling  in  VTB  environ¬ 
ment,  especially  in  the  case  where  a  battery  SFM  model  itself 
is  very  computation  intensive. 


Fig.  8.  Simulation  run  1:  Li-ion  battery  discharge  with  constant  resistance 
load  of  0.3  Q.  See  Fig.  6  for  the  meaning  of  the  legend. 


Fig.  9.  Simulation  run  2:  Li-ion  battery  discharge  with  constant  current  load 
of  2  A.  See  Fig.  6  for  the  meaning  of  the  legend. 

Simulations  of  electrochemical  cycling  processes  using 
battery  models  with  CCVS  interface  are  demonstrated  in 
VTB  environment  (see  Fig.  10).  In  the  simulation,  the  bat¬ 
teries  are  charged  and  discharged  through  a  current  source 
which  is  controlled  by  a  pulse  signal  generator. 

The  simulated  cell  voltage  responses  to  the  pulse  cur¬ 
rent  are  presented  in  Figs.  11  and  12.  As  observed  and  ex¬ 
pected,  the  cell  voltage  decreases  in  discharge  and  increases 
in  charge.  However,  the  first  principle  Li-ion  battery  model 
more  accurately  reflects  the  cell  voltage  responses  than  the 
simplified  NiMH  battery  model  does  in  the  case  of  the  relax¬ 
ation  phenomenon  [13]  which  is  shown  in  Fig.  12  right  after 
the  short  term  discharge  and  charge  process  as  cell  voltage 
slowly  recovers  even  though  no  current  flows  in  the  circuit 
during  those  periods.  The  relaxation  phenomenon  is  caused 
by  the  existence  of  the  concentration  gradient  of  lithium  ions 


Fig.  10.  Schematic  of  battery  simulation. 
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Fig.  11.  Simulated  NiMH  cell  voltage  response  to  the  pulse  current. 


x 


Fig.  14.  Li-ion  battery  electrolyte  concentration  profiles  in  charge  and  relax¬ 
ation  process  of  the  Li-ion  battery  model.  Dashed  line  divides  the  separator 
and  composite  cathode. 


in  the  electrode  and  electrolyte  formed  in  the  discharge  pro¬ 
cess  and  charge  process.  Lithium  ion  concentration  profiles  in 
the  battery  are  presented  in  Figs.  13  and  14  for  the  first  current 
pulse.  Driven  by  the  concentration  gradient,  solution  phase 
lithium  ions  at  the  electrolyte-electrode  interface  which  are 
consumed  or  regenerated  by  the  electrochemical  reactions  are 
replenished  or  dispersed  through  diffusion.  Because  no  diffu¬ 
sion  phenomenon  is  considered  in  the  NiMH  battery  model, 
the  relaxation  behavior  does  not  show  up  in  the  results  as 
expected. 


4.  Conclusions 


Fig.  12.  Simulated  Li-ion  cell  voltage  response  to  the  pulse  current. 


Fig.  13.  Li-ion  battery  electrolyte  concentration  profiles  in  discharge  and 
relaxation  process  of  the  Li-ion  battery  model.  Dashed  line  divides  the  sep¬ 
arator  and  composite  cathode. 


The  resistive  companion  modeling  of  batteries  with  two 
electric  type  terminals  is  introduced  in  the  paper.  The  methods 
can  be  generally  applied  to  other  physical  type  models  such  as 
mechanical  or  hydraulic  ones  so  long  as  energy  conservation 
equations  are  satisfied  in  the  form  of  across  and  through  vari¬ 
ables.  Native  RC  modeling  method  and  the  one  using  CCVS 
interface  are  presented  in  accordance  to  the  complexity  and 
computation  requirement  of  battery  models.  It  is  conclude 
from  our  work  that  using  CCVS  interface  is  currently  the 
first  choice  to  model  batteries  in  VTB.  Simulations  using 
RC  models  with  CCVS  interface  are  demonstrated  on  VTB 
platform  and  some  results  are  presented  and  discussed.  The 
ready  availability  of  the  VTB  software  and  models  of  electric 
equipments  make  the  VTB  environment  useful  to  study  the 
design  and  performance  of  advanced  power  sources. 
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